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o : 

(N , 

In dimension d > 3, the directed polymer in a random medium undergoes a phase transition 

, between a free phase at high temperature and a low temperature disorder dominated phase. For 

. the latter phase, Fisher and Huse have proposed a droplet theory based on the scaling of the free 

' energy fluctuations AF(/) ~ at scale I. On the other hand, in related growth models belonging to 

I the KPZ universality class, Forrest and Tang have found that the height-height correlation function 

, is logarithmic at the transition. For the directed polymer model at criticality, this translates into 

' logarithmic free energy fluctuations AFt^{1) ~ (In/)"^ with a = 1/2. In this paper, we propose a 

I I ' droplet scaling analysis exactly at criticality based on this logarithmic scaling. Our main conclusion 

^ , is that the typical correlation length £,{T) of the low temperature phase, diverges as ln5(r) ~ 

^ • (-ln(Te-T))^/'' ~ (-ln(Tc-T))^, instead of the usual power-law ^(r) ~ {T^-Ty. Furthermore, 

C/3 ' the logarithmic dependence of AFt^{1) leads to the conclusion that the critical temperature Tc 

, actually coincides with the explicit upper bound T2 derived by Derrida and coworkers, where T2 

• ' corresponds to the temperature below which the ratio diverges exponentially in L. Finally, 

I since the Fisher-Huse droplet theory was initially introduced for the spin-glass phase, we briefly 

a mention the similarities and differences with the directed polymer model. If one speculates that the 

I ' free energy of droplet excitations for spin-glasses is also loga rithmic at Tc, one obtains a logarithmic 

1-^ I decay for the mean square correlation function at criticality C^(r) ~ 1/ (Inr)"^, instead of the usual 

pH , power- law l/r'*~^+''. 

o ■ 

' I. INTRODUCTION 

(N : 

^ ■ The model of a directed polymer in a random medium plays the role of a 'baby spin glass' model in the field of 
' disordered systems P, 0, IIHLIM; At low temperature, there exists a disorder dominated phase, where the order 
■ parameter is an 'overlap' 0- In finite dimensions, a scaling droplet theory was proposed in direct 

I correspondence with the droplet theory of spin-glasses whereas in the mean-field version of the model on the 
. Cayley, a freezing transition very similar to the one occurring in the Random Energy Model was found 0I- The 
\^ phase diagram as a function of space dimension d is the following In dimension d < 2, there is no free phase, 
, i.e. any initial disorder drives the polymer into the strong disorder phase, whereas for d > 2, there exists a phase 
transition between the low temperature disorder dominated phase and a free phase at high temperature [l^ 
J2 [ where the free energy has its annealed value. This phase transition has been studied exactly on a Cayley tree 
on hierarchical lattice • In finite dimensions, bounds on the critical temperature Tc have been derived [ill ll^ 
Toid) <Tc< 72 (d). The upper bound 12 (d) corresponds to the temperature above which the ratio Z\l{Zi,)'^ remains 
finite as L — > 00. The lower bound Tq corresponds to the temperature below which the annealed entropy becomes 
Q negative. In c? = 3, the critical properties have been studied numerically bv Il3lll5l |. with different conclusions. The 
O ' study of 13] gives a slightly negative value a ~ —0.1, whereas the work of [l^ yields a correlation length exponent 
~ 4, corresponding through hyperscaling, to a = 2 — v ^ —2. 
. In related growth models belonging to the KPZ universality class, numerical studies and theoretical arguments 
rN ' [1^ IrR have found that the height-height correlation function is logarithmic at the transition. For the directed 
polymer model at criticality, this translates int o log arithmic free energy fluctuations A_F(/) ~ (InZ)'^ with an exponent 
cr = 1/2 that has been measured in d = 3 |l.5l llfi| ■ In this paper, we make a droplet analysis at criticality based on 
this logarithmic scaling, in direct correspondence with the Fisher-Huse droplet theory of the low temperature phase 
based on the free energy scaling AF(Z) ~ . The matching between the two droplet distributions below and at 
allows us to derive that the typical correlation length ^ (T) of the low temperature phase does not follow a power-law 
(Tc - T)-", but diverges instead as In^(T) - (In l/(Tc - Tyfl" . Moreover, we argue that the logarithmic fluctuations 
of the free energy at criticality leads to the conclusion that the critical temperature actually coincides with the 
upper bound Ti derived by Derrida and coworkers, since Ti corresponds to the temperature below which the ratio 
Z£/(Zl)-^ diverges exponentially in L. 

The paper is organized as follows. In Section m we explain on the pure wetting and Poland-Scheraga model of 
DNA denaturation how the transition can be analyzed in terms of the distribution of large loops. The same approach 
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will then be adapted in the remainder of the paper to analyse the freezing transition of the directed polymer, using 
the loop distribution between two independent copies of the polymer in the same disordered sample. In section llTTI we 
describe the Fisher-Huse droplet theory of the low temperature phase based on the scaling AF{1) ^ . In section Hvl 
we describe the droplet theory based on the logarithmic scaling of the free energy at criticality AF(Z) ~ (ln^)°' with 
a = 1/2, and we obtain the divergence of the correlation length £,{T) near Tc, as well as the behavior of the overlap. 
In section Ivl we argue that the transition temperature Tc coincides with the upper bound T2 derived by Derrida and 
coworkers, and discuss why it is not the case in other disordered systems. Finally, in section IVTl we discuss the case 
of spin-glasses for which the Fisher-Huse droplet theory was initially developed. If one assumes, by analogy with the 
directed polymer model, that the free energy of droplet excitations is logarithmic at Tc, one obtains some unusual 
behavior for the correlation function at criticality. Our conclusions are summarized in Section I Vlll In the Appendix 
IXl we explain in more detail the matching procedure for correlation function that we use in this paper. 



II. PURE DELOCALIZATION TRANSITIONS : ANALYSIS IN TERMS OF THE LOOP 

DISTRIBUTION 



A. Wetting and Poland-Scheraga model of DNA denaturation 



The wetting model is defined by the partition function 

exp /3 ^ 

(1) 

RW \ l<a<N J 

where the sum is over one-dimensional random walks (RW) of 2L steps, starting at z(0) — 0, with increments 
z{a -|- 1) — z{a) = ±1. The random walk is constrained to remain in the upper half plane 2 > 0, but gains an 
adsorption energy eo if z{a) = 0. 

The Poland-Scheraga (PS) model of DNA denaturation is closely related to the wetting model. It describes 
the configuration of the two complementary strands as a sequence of bound segments and open loops. Each loop of 
length I has a polymeric entropic weight 

M{1)^^ (2) 

whereas each contact has a Boltzmann weight e~^'^°. The partial partition function Zps(l,a) with bound ends at 
monomers (1) and {a) satisfies the simple recursion relation 

Zpsil, a) = e-^'- -^(^ - «')^P5(1, «') (3) 

a' = l 

The wetting model corresponds to a Poland-Scheraga model with parameter = 2 and loop exponent c = 3/2 
(this exponent comes from the first return distribution of a one dimensional random walk). These models without 
disorder are simple to solve explicitly : they undergo a phase transition between a localized phase at low temperature, 
characterized by an extensive number of contacts, and a delocalized phase at high temperature. The transition is first 
order for c > 2, and second order for 1 < c < 2 with correlation length exponent v — l/(c — 1). Let us now describe 
how this transition can be understood from the point of view of the loop distribution. 

B. Loop distribution in the low temperature phase 

For the wetting or Poland-Scheraga model with loop exponent c, the loop distribution normalized to unity 



jdlP{l,T)^l (4) 



can be explicitly computed |2di . Near Tc, it is useful to decompose it into two terms 

P{UT)^ Pf^rnte il,T)+ Piarge (^, T) (5) 
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The first term represents the statistics of finite loops / = 1,2, .., whereas the second contribution concerning large 
loops I ^ 1 follows the scaling form 



dlPlar,e{l,T)=M{T)j<i>(^-^^] (6) 



Here £,{T) is the correlation length that diverges at the transition ^(Tc) = oo, and the factor Af{T) represents the 
'normalization' of large loops. In Poland-Scheraga model, the dependence in I of the probability of a large loop of 
length / involves the entropic weight 1//^ of a free loop of length / entering the definition of the model (O and the 
extensive free energy cost AF{1) — (where we have used hyperscaling) 



so that the scaling function $ in (0 reads 



Plarge{l,T) ^ ^^e-''^ (7) 



Let us now consider the number of contacts nL(T). In the low temperature phase, it is extensive and simply reads 

= <5>. + </>. 

where < b >t represents the averaged length of sequences of consecutive bound monomers (which remains finite as 
T Tc), and where < I >t represents the averaged loop length of the full distribution |SJ) 

<1>T^ J dll[Pf„,ue{l,T) + Piarae{l,T)]^ finite + j dm{T)^l^-^^ ^ finite + M{T)i{T) j d\^{\) (10) 

For 1 < c < 2, this averaged loop length diverges as 

<1>T-U{T)^{T) (11) 

In wetting and Poland-Scheraga models, the energy is directly proportional to the number of contacts so that the 
energy density also vanishes as 

""^^^ ^ ^ mWi ^^^^ 

On the other hand, the singularity of the energy is given by the derivative of the free energy density f{T) ~ 1/^(T) 
with respect to temperature. The critical behavior of the energy is governed by 

The comparison between Eqs (I12f) and \i'6\i yields the following differential equation for the correlation length £,{T) 
in terms of the normalization J\f{T) 

d'^^i^) = -±- (14) 

C. Loop distribution at criticality 

At criticality, the loop distribution is simply given by the entropic weight that enters the definition of the model 
@ 

dlPlargeil,T,) (15) 



4 



For c > 2, the averaged length < I >= J dl lPiarge{l,Tc) is finite, so that the number n(Tc) of contacts with the 
substrate is stih extensive n{Tc) ^ L and the transition is first order. For 1 < c < 2, the averaged length < I > 
diverges 

< I >T,= J dllPiargeil,T,) = OO (16) 

The number of contacts is not extensive anymore at criticality nL{Tc)/L ^ 1/ < I >Tc= 0. Since the Levy sum of 
n independent variables U drawn from the distribution H15|l scales as Zi + ... + In ^ rt^/^'^^^-', the number of contacts 
^pure J- ^ sample of length L scales as 

nT"(T,) ^ L^-^ (17) 
and the transition is second order. Let us now make the connection with the low temperature phase. 

D. Matching the loop distribution in the critical region 

On length of order I ~ ^ [T) , the two expressions (jSJ and H15(l for the loop distribution for T < Tc and for T = Tc 
should be of the same order 

Plarge{aT),T) = Pfarge (Cl^) , Tc) (18) 

This matching determines the normalization factor JV{T) in terms of the correlation length ^(T) 

Af{T) ^ iUT))'-' (19) 

So as the transition is approached T Tc, the correlation length diverges £,{T) +00, but the density of these large 
loops vanishes N{T) 0. 

The differential equation (|14|l becomes a closed equation for ^(T) 

ffl ^ (aT)r (20) 
The integration with the condition £,{Tc) — oo gives 

f'-\-00 T 

Tc-T^ ^-(e(T))^-^ (21) 



leading to 



C(T) ~ {Tc - T)-"^ with V = (22) 



in agreement with the exact solution. Note that the normalization of large loops vanishes linearly l|19|) independently 
of the loop exponent c 

N{T)^{i{T) f-' ^{Tc-T) (23) 

E. Important ideas for the following sections on disordered systems 

In this section, we have explained how the delocalization transition for the pure wetting or Poland-Scheraga models 
could be interpreted as a vanishing density N{T Tc) — > of large loops of characteristic size that diverges 
^(T — > Tc) — > 00. An important point for the following sections is that the notion of loop distribution which is 
natural to consider in the low temperature phase where the number of contacts is extensive nL{T) L, has still a 
meaning exactly at Tc even if the number of contacts is not extensive anymore. It looses its meaning only in the 
high temperature phase where the number of contacts remains finite. The main idea of the following sections in thus 
that in disordered systems presenting a low temperature phase where the order parameter is an overlap, the droplet 
distribution is not only useful for T < Tc but has still a meaning exactly at Tc even if the overlap vanishes. Moreover, 
it is precisely this critical droplet distribution that determines the properties of the transition. 
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III. DIRECTED POLYMER : LOOP DISTRIBUTION IN THE LOW TEMPERATURE PHASE 

The directed model is defined by the partition function over d dimensional random walks x{i) of L steps 

L 

Zl{I3)^ g-/3£;({:.(i}) ^ith E{{x{i}) = Y€{i,x{i)) (24) 

{x{€)} ^=l 

where the random energies e(i,x) are independent. We have already described in the Introduction its main features 
with the corresponding references. 

A. Statistics of excitations above the ground state 

The droplet theory for directed polymers 0, j is similar to the droplet theory of spin-glasses @ . It is a scaling 
theory that can be summarized as follows. At very low temperature T — > 0, all observables are governed by the 
statistics of low energy excitations above the ground state. An excitation of large length I costs a random energy 

A£:(/) - fu (25) 

where u is a positive random variable distributed with some law Qq{u) having some finite density at the origin 
Qo{u = 0) > 0. The exponent 9 is the exponent governing the fluctuation of the energy of the ground state is exactly 
known in one-dimension 0{d = 1) = 1/3 21, 22, 23, 24] and for the mean-fleld version on the Cayle y tr ee 9(d = oo) = 
0. In finite dimensions d = 2,3,4,5, the exponent 6{d) has been numerically measured [25L l26l [271 l28j . and we 
quote the results of the most precise study we are aware of (2^ for dimensions d = 2, 3 : 9{d = 2) = 0.244 and 
6{d = 3) = 0.186. Note that the existence of a finite upper critical dimension dc where the exponent would vanish 
9{dc) = has remained a very controversial issue between the numerical studies [25II26L l2^ l2^ and various theoretical 
approaches p^l30ll3H. 

From H25|l . the probability distribution of large excitations I ^ 1 reads within the droplet theory 

dlp{E = 0, ^ ^^-PAE(i) ^ dl^-pi^^ (26) 

where the factor dl/l comes from the notion of independent excitations 0. In particular, its average over the disorder 
follows the power-law 

dlp{E^O,l)^ duQo{u)je-^'"-^TQiO)j^ (27) 

This prediction describes very well the numerical data in the regime 1 ^ Z <C i in dimensions d — 1, 2, 3 [3^ . 



B. Low temperature phase governed by a zero-temperature fixed point 

According to the droplet theory, the whole low temperature phase < T < Tc is governed by a zero-temperature 
fixed point. However, many subtleties arise because the temperature is actually 'dangerously irrelevant'. The main 
conclusions of the droplet analysis @ can be summarized as follows. The scaling governs the free energy cost of 
an excitation of length I, provided one introduces a correlation length ^(T) to rescale the length I 

AF{l)=(^-^^\ (28) 

Here as before, u denotes a positive random variable distributed with some law Q{u) having some finite density at 
the origin Q(u = 0) > 0. Moreover, this droplet free energy is a near cancellation of energy and entropy contributions 
that scale as 

AE{1) ~ l^/'^w (29) 

where w is a random variable of order 0(1) and of zero mean. The argument is that the energy and entropy are 
dominated by small scale contributions of random sign j^] j whereas the free energy is optimized on the coarse-grained 
scale (.{T). These predictions for the energy and entropy have been numerically checked in [5ll33|. 
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C. Loop distribution for two polymers in the same disordered sample 



We now describe how the analysis of Section |n] for the pure transition of the wetting or Poland-Scheraga models 
can be adapted to the present disordered case. For T < T^., the number of contacts of two independent polymers x(i) 
and y{i) in the same disordered sample 



riL 



(T) = ^ < > (30) 



is extensive, and the density of contacts, also called the overlap, is precisely the order parameter of the low temperature 
phase 0,11 HQ 

,(T) . hm ( ^) (31) 



L- 



^00 



Note that on the Cayley tree where = 0, the distribution of this overlap is made of two delta peaks at g = and 
at g = 1 0, whereas in finite dimensions with 9 > 0, the distribution of this overlap is expected to be a single delta 
fimction at q{T) 

One may thus analyse the configuration of two polymers in the same sample in terms of contacts separated by 
loops. Again, it is useful to make the decomposition ((Sj to concentrate on the contribution of large loops that follows 
a scaling form based on the free energy scaling of a droplet of length / (|28|) 

dlPiargeil^ T) = AA(T) | e'^^^^') = ^(^T)je-^'iimT ^ (32) 

As in the wetting case, the factor Af{T) represents the normalization of large loops that will be determined below. 

The important difference with respect to the pure case is that now the probability of a large loop at a given spatial 
position depends upon a random variable u. This introduces very different behaviors for typical and averaged loop 
distributions. The averaged loop distribution has the following power-law decay in the whole low temperature phase 



dlPlargeil) ^ ^MiT) [ duQ (u) ^ M(^T)Q{0)T^ ( ^) (33) 



whereas the typical decay is an exponential with exponent 9 



In Piarge{l)^HmT)/ 1) ~f3(^^] Uo (34) 



with uo = Jq°° duuQ{u). In particular, the typical distribution 



dl 

diF''"''{L) =M(;r; 

has a finite first moment 



dZP*''f(/) =AA(T)ye"''(«rF))''"'' (35) 



< I >large^ / dllP'yP{l) ^ N{T)i{T) (36) 



As the transition is approached, this term will govern the divergence of the first moment of the full droplet distribution 

< / >tot^< I > f.mte + < I >large- AA(T)e(r) (37) 

As a consequence, the number nL{T) of contacts (and equivalently the number of loops) will have a vanishing density 
of order 

L <1>T AA(r)C(T) ^ ' 
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D. Energy fluctuations near the transition 

Let us now consider the specific heat c(T) that measures the fluctuations of the energy 

T^c{T) = -< {El- < El > f >= ^ ' 



(39) 



In the high temperature phase, the specific heat coincides with its annealed value. In particular, for the model H24() 
with the following Gaussian distribution for the random energies e{i,x) 



the annealed specific heat reads 



Pie) 



^ ^anniT^ — 1 



(40) 



(41) 



meaning that the energy fluctuations per monomer are Ae^ = 1. In the low temperature phase, the free energy 
is above the annealed free energy, and this implies that the specific heat is less than the annealed value, with a 
non-diverging singular part 3j. More precisely if one introduces the specific heat exponent a, we have 



c{T) - Cann{Tc) A{Tc - T)-"" with A>0 and a<0 



(42) 



To interpret this loss of energy fluctuations in the low temperature phase in terms of the droplets, it is convenient 
to write the energy fluctuations in terms of the energy difference between two copies x{i) and y{i) of the polymer in 
the same disordered sample 

1 < [El- < El >)' >= ^ < (E (M^)}) - E ({y(*)}))' ^ < ^(*)) > (^3) 



i=l 



l<i<j<L 



In this sum, the monomers i corresponding to a contact x{i) = y{i) do not contribute, whereas the monomers i inside 
loops x{i) 7^ y{i) are expected to have an energy fluctuation of order Aef ^ 1 fEa 1^ . As explained above, for 
T > Tc all monomers are characterized by energy fluctuations Aef ~ 1. We thus write that the extensive loss in 
energy fluctuations just below Tc is proportional to the density of contacts x{i) = y{i) obtained in Eq (|38|l 



< {El- < El >f > 



< {El- < El >f > 



L 



n{T)^{T) 



Near T^, the free energy density /(T) = —Th\ZL{0)/L behaves from hyperscaling as 

1 



/(T) - /(T,) = 



Deriving twice with respect to temperature H39() . the singular part of the energy fluctuations behave as 



- < {El- < El >f > 



- < {El- < El > f > 



(f f 1 



dT2 \^{T) 



(45) 



(46) 



(47) 



The consistency of Eqs H45|l and (j47|l then gives the following differential equation for the correlation length ^ (T) 
in terms of the normalization M{T) 



f 1 



dT2 Ve(r) J Af{T)C{T) 



(48) 
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IV. DIRECTED POLYMER : LOOP DISTRIBUTION EXACTLY AT CRITICALITY 
A. Logarithmic fluctuations of the free energy 

Let us now consider what happens for T = Tc- Forrest and Tang fl^ have conjectured from their numerical results 
on a growth model in the KPZ universality class and from the exact solution of another model 34] that the fluctuations 
of the height of the interface were logarithmic at criticality. For the directed polymer model, this translates into a 
logarithmic behavior of the free energy fluctuations at 

AF{L,Tc)^ (InLYv (49) 

where ?; is a positive random variable of order one distributed with some law R{v), and where the exponent was 
measured to be in d = 3 0, 0| 

a = i (50) 

Further theoretical arguments in favor of this logarithmic behavior can be found in Il8| . The argument of 
is that the power-law behavior F(L, Tc) ~ L^" is impossible at criticality so that 9c = 0. From the scaling 
relation 6c = 2C,c — 1 between exponents [s^ , the roughness exponent C is expected to be exactly Cc = 1/2 HjI , and a 
renormalization argument then leads to logarithmic fluctuations of the free energy . 

B. Loop distribution at criticaHty 

With the scaling H49|) . the loop distribution exactly at criticality reads 

dlP^^{l) = ^e-/3AF(0 ^ rfig-/3.(lnO-. (51) 

The averaged loop distribution has then the following extremely slow decay 

i Jo i(lnt)'^ 

whereas the typical decay is given by 



lnPTe(0 = -ln^-/3c(lnO''wo (53) 
with vo = J^°° dvvR{v). The typical distribution 

dlP^lP{l) = ^6-^(1"')'"° (54) 
has an infinite first moment for < tr < 1 (independently of the prefactor /3c wo) 

<1>T^= J dllP^yf = oo (55) 
and the number of contacts ulITc) is not extensive in L. 

C. Matching the typical loop distribution in the critical region 

For I ~ ^, the two expressions of the typical correlations in the low temperature phase and at criticality 
should coincide as in (|18() 

Ptyp{aT),T)^Ptyp{aT),Tc) (56) 
This gives the following relation between the normalization and the correlation length ^(T) 

AA(T) ~ e-^P'^«)" (57) 
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where K ~ Pc^o is some constant. 

Using (|48l) . we thus obtain the following closed differential equation for the correlation length £,{T) 

Near the boundary condition ^(Tc) ~ oo, the leading divergence of ^(T) is given by 

1 ,K(lnO 



e^^''^?) (59) 



We thus obtain at leading order that the correlation length ^(T) diverges as 

e(T) - e(*'"T^)'^"+- (60) 
instead of the usual power-law behavior ^{T) ^ (Tc — T)^^ . The free energy difference thus vanishes very slowly as 

/(T) - /(T,) - ^ - e-(*'"^)'''+- (61) 

The normalization of large droplet vanishes at leading order as ()48|l 

7V(T) - {Tc - Tf (62) 

and the contact density or overlap as (|31|1 



In particular, for the value a — 1/2 measured in d = 3 [15|,ll6j, we obtain the final results 

f{T)-HTc) (64) 
q{T) ~ e-(*'"'J^)'+2>"'Jv^+- (65) 

V. DIRECTED POLYMER : EXPLICIT VALUE OF THE CRITICAL TEMPERATURE 

A. Exact bounds on Tc derived by Derrida and coworkers 

Let us first recall the physical meaning of the exact bounds for the critical temperature derived by Derrida and 
coworkers [Tl lll^lll | 

To(d) < < T2{d) (66) 
The upper bound T2{d) corresponds to the temperature above which the ratio 

Rl{T) = (67) 

remains finite as L ^ oo. The lower bound Tq corresponds to the temperature below which the annealed entropy 
becomes negative. 

In dimensions d = 1,2, the upper bound is at infinity T2 — 00, whereas for c? > 3, the upper bound T2 is finite. 
The interpretation is as follows (1J|]. The ratio H67I) can be decomposed according to the probability PL{m) that two 
independent usual random walks in dimension d meet m times before time L 

L 

Rl{T) = Pi"^)B"' (68) 
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where the factor 

can be expHcitly computed for any distribution of the site disorder variable e. In dimensions d — 1,2, two random 
walks meet an infinite number of times as L ^ oo, whereas for d > 3, they meet a finite number m of times as L — > oo. 
The distribution of m decays exponentially 

P(m) - (1 - A)A"' (70) 
where (1 — A) is the finite probability of never meeting again. T2 is defined as the temperature where 

AB{T2) - 1 (71) 
For T > T2, BiT) < B{T2) = l/A, and the ratio Rl{T) has a finite limit 

Roo(T>T2) = — ^7"^ ^ (72) 

For T <T2, Rl{T) is a geometric series of parameter AB{T) > 1, and it thus diverges exponentially in L 

L 

Rl{T < T2) - (1 - A) ^ (AB(T))'" - {AB{T))^ (73) 

m— 1 

Exactly at T2, the ratio diverges but not exponentially 

L 

Rl{T2) = {1-A)Y,I^L (74) 



B. Interpretation in terms of the probability distribution of free energies 

Let us now interpret the above results of the ratio i?^ (T) in terms of the probability distribution {F) of the free 
energy F — — kT In Zl over the samples of length L. By definition H(i7l) . one has 

, , J dFPL{F)e-^'^^^ . , 

iJdFPUF)e-0^-y (^^^ 



For T > T2 (d) , the ratio Roc (T) is finite : this means that the fluctuations of the free energy over the samples 



i^PL]Lmpies = I dFPUF)F' - (J dFP^{F)F 



(76) 



remain of order 0(1) in the limit L — > cxd. 

On the other hand, for the directed polymer in the low temperature phase T < Tc, the fluctuations of free energies 
over the samples is expected to have the same scaling as the fluctuations of free energies within the same sample when 
the end-point varies |5| : the fluctuations of free energy over the samples are thus governed by the droplet exponent d 

[Ai^L],_,,e, (T < Te) ^ [APi]„,„^,^, (T < T,) ^ L'^ {11) 

Let us now recall Zhang's argument Q that allows to determine the exponent 77 of the tail of the free energy distribution 

Pl{F ^ -00) ^ e^'^iy) (78) 

Moments of the partition function can be then evaluated by the saddle-point method, with a saddle value F* lying 
in the negative tail 178() 

Zf = J dFPLiF)e-'"^^^ ~ J dFe-^T^y e-"/^^^ ^ (79) 

Since these moments of the partition function have to diverge exponentially in L, the exponent rj of the tail H78|) reads 
in terms of the droplet exponent 

(80) 
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C. Why Tc coincides with T2 for the directed polymer in finite dimensions 

At Tc, the fluctuations of the free energy are expected to be logarithmic, as discussed around Eq 

AF-(hii)'" with cr==i (81) 

From these logarithmic fluctuations, it seems rather difficult to obtain an exponential divergence in L of the ratio 
Rl{Tc) (|75() . so the strict inequality < T2 seems very unlikely. On the contrary, if Tc = T2, it is very natural to 
obtain the divergence found for the ratio at T2 (fZ^ 

i?L(T2)^T~e'"^ (82) 

Moreover, to obtain the linear divergence (|74|l . the saddle-point method described above for the low temperature 
phase H79f) gives that the tail of the free energy distribution should be at criticality 

-PTe(T -00) - e~(o^) " with T^c = (83) 

1 — cr 

The value (7 = 1/2 corresponds to the tail exponent rjc = 2. 

Our conclusion is thus that the critical temperature in finite dimension d coincides with the temperature T2{d) 

T,{d)^T2{d) (84) 

Explicit expressions for T2{d) in terms of usual integrals appearing in the theory of random walks can be found in 
[III [13 for site and bond disorder respectively. 

As a final remark, let us mention that for the directed polymer on the Cayley tree that plays the role of a mean-field 
version of the model, the critical temperature does not coincide with the upper bound T2, but coincides with the lower 
bound To H()()|l below which the annealed entropy becomes negative |lll] 

Tc{Cayley) = To{Cayley) < T2{Cayley) (85) 

This shows that the mean-field limit of the tree structure changes the nature of the transition with respect to the 
finite-dimensional case. The technical reason seems to be that 6' = and AT — 0{\) in the low temperature phase 
of the Cayley tree, whereas Zhang's argument above is consistent only if r; = 1/(1 — ^?) > 1 to ensure the convergence 
in the presence of the exponential term e~^^^ (F,a I79|l . When 9 — 0, the tail of Pl{F —00) is also an exponential 
^aF/{AF) jj-^ ^Yie Random Energy Model and one has to take into account the minimal free energy that can be 
obtained for a finite size L. From a physical point of view, the reason is that the configurations of two polymers in 
the same disordered sample are very different. In finite dimensions, contacts and loops alternate extensively, whereas 
on the tree, the loops simply do not exist : the two polymers may only coincide over some distance and then never 
meet again. Since the exponential tail found for the Cayley tree actually corresponds to the universal Gumbel tail 
for the minimum of independent variables, this shows that the non-exponential tail found in finite dimensions reflects 
the importance of correlations between the free energies of paths due to the presence of loops. 



D. Why Tc is different from T2 in other finite dimensional disordered systems 

The fact that the fluctuations of free energies over the samples have the same scaling as the droplet excitations 
within one given sample (|77|l is very specific to the directed polymer model. In other disordered models, such as 
spin-glasses for instance, the fluctuations of free energies over the samples scale instead as js^ [s^l . 

at any temperature. This scaling simply reflects the Central-Limit fluctuations of the L'^ disorder variables defining 
the sample. The directed polymer escapes from these normal fluctuations because it is a one-dimensional path living 
in a 1 + d disordered sample : each configuration of the polymer only sees L random variables among the T^+'' disorder 
variables that define the sample, and the polymer can 'choose' the random variables it sees. So for other disordered 
systems having fluctuations over the samples governed by (|86ll . the ratio Rl{T) will diverge exponentially as any 
temperature. The temperature T2 is thus infinite 



T2 = 00 



(87) 
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and has nothing to do with any critical temperature. However, the droplet exponent 9 is expected to govern the 
correction to the extensive part of the mean value 

Fl ^ L^'fo + LVi (88) 

It can for instance be measured in the free energy difference upon a change of boundary conditions that forces the 
introduction of some domain wall in the sample j9|. 

VI. SPIN-GLASSES 

In this Section, we briefly mention the analogies and differences between spin-glasses and the directed polymer 
model described above. 

A. Fisher-Huse droplet theory of the low temperature phase 



Discussions on the droplets statistics in pure Ising models can be found in Refs |38ll83 . l4r| . Here we only summarize 
very briefly the Fisher- Huse droplet theory of the spin-glass phase Q . 
The free energy to make an excitation of large size I follows the scaling 



AFil,T)^ i^-^j u (89) 

where m is a positive random variable of distribution Q{u). 
The truncated correlation function 

C(r) - I < SoSr > - < So >< Sr > \ (90) 
for large distance r is governed by the probability that the two points o and r belongs to the same large droplet 

C(r) -AA(r)e"(«(^)'" (91) 
where the prefactor JV{T) corresponds to the Edwards- Anderson parameter 

AA(r) ~ am) - QEAiT) = ^ E < >' (92) 

i 

The typical decay of the correlation function is 



lnC{r)^\nqEAiT)-{ ^ ) uo (93) 



whereas the mean-square correlation function 



C^{r)^TQ{0)qlAT)[^^j (94) 
has a power-law decay in the whole low temperature phase. 

B. spin-glasses at criticality 

As mentioned above, the exponent of excitations is expected to govern the correction to extensivity of the averaged 
free energy in the low temperature phase (Eq. I88|l . At criticality where the extensive part of the free energy vanishes, 
the averaged free energy is thus expected to be governed by the free energy scale of excitations at Tc 

roplet 

(95) 
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According to the argument of [T3 , a power-law behavior 

niT,) ~ L'^ (96) 

is not possible at a critical point, so that 9c — 0- The droplet free energy at criticality can thus only be logarithmic 

[^FLin)],,„p,,,^{lnLrv (97) 

where w is a random variable. In pure models, the power-law decay of truncated correlation can be interpreted as a 
logarithmic cost with apure = 1. It therefore seems natural to expect < cr < 1, since in the low temperature phase, 
the free energy cost in the presence of disorder is much less than the free energy cost l"^"^ of pure systems. 
Then the truncated correlation H90() is expected to behave at large distance as 

C(r) ~ e-'^^^-^^^) ^ e-'^^i"'-)"" (98) 

The typical decay is then given by 



InC(r) =-(lnr)''/3t)o (99) 
whereas the mean-square correlation function decays very slowly as 

TQ(0) 



CHr)^^ (100) 



in contrast with the usual power-law correlation C^(r) ~ i/r'^-^+v g] xhe matching at scale r ~ £,{T) between the 
typical correlation functions at T < (Eq. I93|l and at (Eq. Ref EH ) yields the following relation between the 
Edwards- Anderson parameter qsAiT) and the correlation length ^{T) 

lnqEA{T)^-{lnaT)r (101) 

For the directed polymer model, it was possible to carry the discussion further because the one dimensional structure 
provides a connection between the divergence of large loops and the total number of small loops (|^ . In spin-glasses, 
the relation between properties of large droplets that dominate the correlation function and properties of small droplets 
that dominate the energy fluctuations depends on geometric assumptions on the shape of the droplets and on their 
spatial organization. This goes beyond the present work. 



VII. CONCLUSIONS AND PERSPECTIVES 



In this paper, we have proposed a coherent picture of the freezing transition of directed polymers in dimension 
d > 3 from the following building blocks 

(i) the Fisher-Huse droplet theory of the low temperature phase 

(ii) the Forrest- Tang result concerning the height-height correlation exactly at criticality in related growth models 
belonging to the KPZ universality class 

(iii) the exact bounds on the critical temperature of Derrida and coworkers. 

Our main conclusions are that the critical temperature Tc coincides with the upper bound T2 derived by Derrida 
and coworkers, and that the logarithmic fluctuations at criticality AF{l,Tc) ~ (InZ)"' with a — 1/2 leads to unusual 
critical properties. In particular, the typical correlation length ^(T) of the low temperature phase, diverges as 
InC(T) - (-ln(Tc - T)y/'' - i-\n{Tc ~ T)f, instead of the usual power-law ^(T) - (T^ - T)"". These results 
emerge from the following picture. Below Tc, the number of contacts of two polymers in the same disordered sample 
is extensive, and their configuration can be thus described as a sequence of contacts and loops. As the transition is 
approached, the averaged length loops diverges, but our main point is that the loop distribution has still a meaning 
exactly at Tc. (It looses its meaning only in the high temperature phase). Moreover, it is precisely this critical droplet 
distribution that determines the properties of the transition. 

Finally, since the directed polymer plays the role of a baby spin-glass model, we have briefly mentioned some 
similarities and differences with the directed polymer model. We have discussed some consequences for the Edwards- 
Anderson order parameter and the mean-square averaged correlation, if one speculates that the free energy of droplet 
excitations for spin-glasses is also logarithmic at Tc. Note that the freezing transition of the directed polymer model 
is completely asymmetric w.r.t. Tc'. there is no singularity in thermodynamic quantities as T —^ T^ ,3,], since the 
free energy coincides with its annealed value for T > Tc- For the spin-glass transition, this raises the question of the 
relations between the critical properties below and above Tc- 

We are presently studying numerically various aspects of the freezing transition of the directed polymer in d = 3 [4H 
to see if the scenario proposed in the present paper can be discriminated from the usual power-law critical behaviors 
used previously in the literature to analyse the data 0, . 
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APPENDIX A: MATCHING PROCEDURE FOR CORRELATION FUNCTIONS IN THE CRITICAL 

REGION 

As suggested by the referee, we explain in more details in this Appendix the matching procedure used in the text. 
We illustrate it with explicit examples concerning the truncated correlation function at large distance 

Crir) =< SoSr > ~ < So >< > (Al) 

in pure spin models. 

1. Finite-size scaling theory for the correlation function 

It is useful to introduce the ratio p — r/^{T) to rewrite the correlation function as 

CT{r)^G(^r.p^^^ (A2) 

At criticality where S,(Tc) — oo, the correlation function is given by the limit p 

CTAr)=G{r,p = Q) (A3) 

In the opposite regime where p ^ 1, the correlation is expected to become a scaling function of the ratio p with a 
temperature dependent prefactor J^{T) 

C.M^^^AA(r)$(,= ^) (A4) 

The matching procedure between the two regimes p — > and p 3> 1 consists in requiring that the two expressions 
should have the same order of magnitude at the matching value p* = r* /^(T) ^ 1 : this determines the normalization 
factor 7V(T) as 

UiT) ^ CrAr* - m) (A5) 
We now illustrate the above scheme with Ising and XY models in d = 2. 

2. The Ising model in d = 2 

Exactly at Tc, the correlation function is given by the power-law of exponent 77 = 1/4 

CrAr)^^ (A6) 
For T > Tc, the correlation function follows the Ornstein-Zcrnickc form [40| 

Ct>t. (r) ~AA+(r) (^^) ^ e"«^ (A7) 
whereas for T <Tc the exact decay at large distance is ^2 



r 



2 



CT<TAr)'^^~{T) ( ) e"^«<% (A8) 



The physical interpretation of this form in terms of droplets can be found in [33, l39l |40j . The matching relation yields 

AA± (T) ~ Ct^ (r* ~ ^ (A9) 

Using ^(T) ^ \Tc — T\^^ , the normalizations of ljA7|l and (|A8p are in agreement with the exact results of [4^ . 
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More generally in spin models in dimension d ior T < Tc, the normalization Af (T) is directly related to the order 
parameter m{T) =< Si(T) > via 

Af-{T)^m^{T) (AlO) 
The decay of the correlation function at criticality defines the exponent 77 

Cr.ir)^^^ (All) 



The matching procedure (|A5|) then reads 

AA(r) ^ CrAr* - ar)) - ^r^^ (A12) 

With Eq (|A10|I . this matching condition simply corresponds to the usual relation 2(i = [d — 2 + vi)v between critical 
exponents. 

3. The XY model in d = 2 

Since the correlation length is infinite in the whole low temperature phase (T < Tc), we consider only the matching 
procedure for T > Tc- Exactly at the correlation function is given by 

CrSr)^^-^ (A13) 



For T > Tc, the correlation function follows the Ornstein-Zernicke form [4 



Ct>t. (r) ~AA+(r) (^) ' (A14) 



The matching procedure (|A5|) yields 

N+{T) ~ CrAr* - m) - ^ ^,},^^^ (A15) 
where the correlation length presents the essential singularity divergence ^{T) ~ e^^^'^''~'^'^ ' . 

4. Discussion 

The above examples show that whenever the critical correlation function is a pure power-law {CT^{r) ^ t'^'Oi the 
off-critical correlation function can be rewritten in the factorized form 

Crir) ^ CrAr)^ I^P = (A16) 

that corresponds to the form ((X4|) in the regime p > 1 with ^'(p) = p''"2+''$(p) and J\f{T) = l/ g'^~^+''(T). 

However, if Ct^ (r) is not a pure power-law, as in the XY case , the simple factorization (jA16|) cannot be written. 

The relations with the polymer models discussed in the text are as follows. The loop distribution P{l,T) plays the 
role of a correlation function. 

In the pure Poland-Scheraga model where the loop distribution at criticality is a power-law p5l) . the off-critical 
loop distribution can be written in the factorized form ^ 

Plarge{l,T) ^ Piargeil ^T^)^ (^^) (Al7) 

as in the Ising model above. On this form, the limit ^(T) ~f oo can be taken to recover the critical distribution. 
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In the directed polymer model where the critical loop distribution is not a power-law but involves logarithm H51|l . 
the off-critical loop distribution (|35|l cannot be written in a factorized form involving the critical loop distribution, 
as in the XY model discussed above. In particular, on the explicit expression H;-{5(l with H57|l . valid in the regime 
I ^ ^{T), one cannot take blindly the limit ^(T) ^ oo to recover the critical loop distribution. 
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